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Abstract 

Motivated by the interpretation of the recent results on the TeV gamma radiation from the Galactic 
center, including the new 2004 HESS data, as a by-product of dark matter particles annihilations, we 
address the question of the largest possible neutralino masses and pair annihilation cross sections in 
super symmetric models. Extending the parameter space of minimal models, such as the mSUGRA 
and the mAMSB scenarios, to general soft SUSY breaking Higgs masses gives access to the largest 
possible pair annihilation rates, corresponding to resonantly annihilating neutralinos with maximal 
gaugino-higgsino mixing. Adopting a model-independent approach, we provide analytical and numer- 
ical upper limits for the neutralino pair annihilation cross section. A possible loophole is given by the 
occurrence of non-perturbative electro-weak resonances, a case we also consider here. We then show 
that a thorough inclusion of QCD effects in gluino (co-)annihilations can, in extreme scenarios, make 
neutralinos with masses in the hundreds of TeV range, well beyond the s-wave unitarity bound, viable 
dark matter candidates. Finally, we outline the ranges of neutralino masses and cross sections for 
models thermally producing a WMAP relic abundance, thus providing reference values for "best-case" 
indirect SUSY dark matter detection rates. 



1 Introduction 



The particle physics nature of the non-baryonic and non-luminous component of the matter budget 
of the Universe stands today as one of the greatest puzzles in the understanding of Nature. A few 
features of this elusive constituent are relatively well known, including its cosmological abundance. 
Structure formation favors a Cold Dark Matter (hereafter CDM) particle candidate, i.e. a particle 
which was non-relativistic at the time of its freeze-out from the thermal bath in the early Universe. 
Combining a wealth of recent observational data, including the first year results on cosmic microwave 
background anisotropies from the WMAP satellite, and assuming a flat ACDM cosmology, the CDM 
abundance has been determined to be, at 95% C.L., [1] 

Ocom/i^ = 0.113 ±0.009 (1) 

where h = 0.71 it 0.04 is the Hubble constant in units of 100 km s~^ Mpc~^ . 

An attractive class of particle candidates for CDM is that of weakly interacting massive particles 
(WIMPs) (see e.g. Ref. [2-4] for recent reviews). One of the bonus of the WIMP scenario is the 
appealing idea that Dark Matter can be a thermal leftover from the early Universe, providing a 
possibly natural and fundamental (in the elementary particle physics sense) explanation for Eq. 
On the other hand, WIMPs can be, in principle, detected, either directly, measuring the recoil energy 
from elastic scattering on nuclei (see e.g. [5] for a recent review), or indirectly, looking for the yields 
of WIMP pair annihilations. In this latter case, privileged targets are neutrinos from the center of 
the Earth or of the Sun [6], natural gravitational traps for WIMPs, gamma rays from spots where 
dark matter might have a large density {e.g. the center of our Galaxy) [7-9], antimatter produced 
in the galactic halo [10-12], or comprehensive multiwavelength analyses which take into account all 
stable products of WIMP pair annihilations as well as their secondary yields (see e.g. the recent 
analysis of Ref. [13] for the case of the Coma cluster). The estimate of the fluxes of the mentioned 
products of WIMP pair annihilations depends on the product of the thermally averaged WIMP pair 
annihilation cross section in the non-relativistic limit, which will be hereafter indicated as {crv), times 
(a suitable integral of) the square of the local number density of WIMPs. The latter quantity is 
inversely proportional to the WIMP at a given Dark Matter density, hence the relevant 

quantity for WIMP indirect detection is (av) /m'^^^p. 

Given a fundamental theory which provides a viable WIMP candidate, one of the theorist's primary 
goals is to understand the overall theoretically viable range of a few basic quantities, such as the WIMP 
mass, its scattering cross section off matter and the above-mentioned ratio of the pair annihilation 
cross section over the WIMP mass squared. A desirable, complementary requirement is that the WIMP 
under consideration thermally produces a relic abundance within the range of Eq. (^. The resulting 
theoretical information is, needless to say, crucial, e.g., to the design and prospects forecast of WIMP 
search experiments, or to the interpretation of experimental results in terms of a WIMP-induced signal. 

Among the few well-motivated extensions of the Standard Model (SM) of particle physics encom- 
passing a WIMP candidate, the option of low-energy supersymmetry has attracted an extensive, and 
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privileged, investigation (see Ref. [14] for a review). The minimal, i?-parity conserving, supersymmet- 
ric extension of the SM (MSSM) offers an ideal WIMP candidate, the lightest neutralino. Contrary 
to other frameworks, for instance the DM candidates of universal extra dimensional scenarios, the 
viable realizations of the MSSM offer a wide range of outcomes for the above mentioned DM particle 
properties, reflecting the a priori lack of information on the supersymmetry breaking mechanism, and 
hence on the supersymmetry breaking lagrangian (see Ref. [15] for a review). 

Restricting our analysis to the theoretical laboratory of supersymmetric models featuring a neu- 
tralino X &s the lightest supersymmetric particle (LSP), we apply here the program of systematically 
investigating the resulting maximal mass and pair annihilation cross sections {o'v). As a con- 
crete, worked-out case study, we assess the possibility of explaining the recent atmospheric Cherenkov 
telescopes (ACT) data on the gamma ray flux from the Galactic center in terms of well defined and 
motivated neutralino DM setups, a scenario which involves particularly large values for both and 
(av). 

In greater detail, we analyze in Sec. I^J the top-down interpretation of the Cangaroo-II and of the 
HESS data on the high energy gamma rays flux from the Galactic center, nailing down both the 
statistically preferred WIMP mass range and the WIMP pair annihilation cross section for a given 
central Dark Matter density. We show, in particular, that the new HESS data from the 2004 campaign 
highly restrict the class of WIMP candidates which provide statistically viable fits to observations. 
In Sec. 13.11 we address the question whether minimal, benchmark, gravity and anomaly mediated 
supersymmetry breaking models are suitable to explain the ACT data, and which are the regions, 
on the {niy^, {o'v)) plane, which those models span. We then point out, in Sec. 13.21 that minimal 
extensions of these benchmark models, involving non-universalities in the Higgs soft supersymmetry 
breaking masses, allow to largely extend the maximal neutralino masses compatible with a thermal 
relic population of neutralinos making up the CDM. Further, those minimal extensions feature large 
annihilation cross sections, and, remarkably enough, the heaviest neutralinos will most likely be within 
reach of future direct DM search experiments. The following Sec. ^ is devoted to a detailed model- 
independent numerical and analytical discussion of the maximal neutralino pair annihilation cross 
section in the MSSM. The latter is achieved in correspondence to resonant annihilation processes of 
maximally higgsino-gaugino mixed neutralinos, with the possible caveat of special, model dependent 
values of the neutralino mass, for which resonant non-perturbative electro- weak effects [16] can produce 
even larger cross sections. The resulting gamma rays spectral function, in this case, is however shown 
to be unfit to explain the HESS data. In Sec. [5] we point out that the inclusion of a non-perturbative 
QCD treatment of the gluino pair annihilation cross section can largely affect the relic abundance of a 
co-annihilating neutralino LSP. In certain scenarios, the explicit violation of the s-wave unitarity limit, 
in models with gluino coannihilations, yields the possibility that neutralinos as heavy as hundreds of 
TeV give a thermal relic abundance in the range of Eq. The concluding Sec. |H1 gives an outlook on 
the maximal cross sections and ^''supersymmetric factors" (o"f)/m^ in the MSSM, both for low-relic 
density models and for models with neutralino relic abundances in the WMAP preferred range. A 
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concise summary of the results presented in this note is provided in the final Sec. [Tj 

2 TeV gamma rays from WIMP annihilations: a model-independent 
analysis 

Two ground based ACTs located in the Southern Hemisphere recently reported the detection of TeV 
gamma-rays from the direction of the Galactic center. The Cangaroo-II telescope observed, during 
2001 and 2002, a very soft spectrum, with a spectral index a ~ —4.6 [17]. The HESS collaboration 
measured, during Summer 2003, when two of the four telescopes were operational, a substantially 
different flux from the Galactic center, featuring a spectral index a ~ —2.2, and achieving a remarkably 
better angular resolution [18]. The results of the 2004 HESS campaign of further observations of the 
Galactic center with the complete array of four imaging ACTs were reported this year in various 
conferences [19-21], and essentially confirm, with greatly improved statistics, the 2003 data. Further 
observations with Cangaroo-III will certainly help clarify the situation, which might be ascribed to 
the presence of multiple sources not resolvable by single telescope observations, or which may depend 
on peculiar instrumental issues, maybe the same which caused other reported discrepancies for various 
different steady sources observed by both ACTs [22]. A further reason could be a time variability, on 
a timescale of around one year, of the observed source; an option which appears highly unlikely, since 
none of the experiments detected any significant source variability [18,20]. 

A few "conventional" astrophysics models have been proposed to explain the very high energy 
gamma rays spectrum observed by ACTs in the center of the Galaxy, ranging from physics involving 
the central supermassive black hole [18, 23] to the production of gamma rays from the interaction of 
accelerated protons (possibly injected by the supernova remnant Sgr A East [24]) with the ambient 
matter [18,25]. Another option which has been investigated is that of the pair annihilation of Dark 
Matter particles living in the proximity of the central regions of the Galaxy, giving raise to a gamma-ray 
continuum from the subsequent decays of the particles' final state products. 

In particular, the latter possibility was pursued in Ref. [26] where the Cangaroo-II and the Whipple 
data from the VERITAS collaboration [27] were analyzed in a model-independent approach using 
an analytical approximation [8, 9] to a putative gamma-ray spectrum from a DM particles mainly 
annihilating into gauge bosons. A correlation between the Dark Matter density in the central region 
of the Galaxy and the DM particle pair annihilation cross section was also presented, at given values 
of the DM particle mass. A similar analysis, again based on the same analytical approximation to the 
gamma-ray spectrum generated by DM pair annihilations, was carried out in Ref. [28] for the 2003 
HESS data. Although no particular fully motivated particle physics setups were considered in those 
papers, since minimal supersymmetric models predict neutralino masses typically much lighter than 
the preferred mass range needed to fit the HESS data, as determined in Ref. [28], a novel, special 
particle physics scenario was proposed in Ref. [29]. The latter features a stable messenger state with 
masses in the tens of TeV range as the DM candidate, in the context of an extended version of gauge 
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mediated supersymmetry breaking models involving an additional Higgs singlet. The option of Kaluza 
Klein DM was finally considered in Ref. [30]. The analysis of Ref. [30] also included contributions 
from internal bremsstrahlung. Although the lightest Kaluza-Klein particle (LKP) fails to produce 
a sufficiently low thermal relic abundance and a large enough LKP pair annihilation cross section 
(without invoking huge boost factors) at particle masses relevant to fit the full gamma rays energy 
range spanned by the 2003 HESS data, it was shown that an hypothetical particle with unsuppressed 
couplings to charged leptons and with increased gauge couplings might give a fairly good fit to the 
data. In this scenario, however, a central Dark Matter density at least a factor 1000 larger than 
what predicted by the cuspy Navarro Frank and White profile [31] would still be needed. Correlations 
between the ACT data and the Egret data [32], and the possibility of a combined DM annihilation 
interpretation of the 2003 HESS data and of the Cangaroo-II data were also addressed in Ref. [33,34]. 
No specific particle physics models were however proposed there to account for the large masses and 
annihilation cross sections to be invoked to explain the HESS data. In Ref. [21] the 2004 HESS data 
were analyzed in terms of supersymmetric or Kaluza-Klein (KK) DM annihilations. None of these 
models was found to give satisfactory fits to the 2004 data, and model-dependent limits were set on 
the WIMPs annihilation cross sections under the assumption of a best-fit power-law background. 

In full generality, the continuum gamma rays spectrum generated by the pair annihilation of a 
WIMP X can be cast as [26] 

d$^ _ „ 19 _9 _i / . „s dA''/ \ , / TTlv \ ~2 



where the symbol / refers to any final state of the WIMP pair annihilation process, yielding a gamma 
rays spectral function (differential number of photons per WIMP annihilation) dN^^ /dE^, and where 
we introduced the quantity (crJ), defined as 

(C^J) ^ ( o^j-2j 3-i ) X ^^(^' (3) 

V 3 X 10 ^''cm'^s ^ / 



where J(Ai}) is the average, over the solid angle AJl, of the following line of sight (l.o.s.) integral over 
the Dark Matter density 

In the following analysis, we assume a solid angle Aft ~ 5 x 10~^ for the Cangaroo-II telescope, 
corresponding to an angular resolution around 0.2", while we consider, for HESS, an angular resolution 
of 5.8', corresponding to the instrumental 50% containment radius^ [18]. The latter superior angular 
resolution will give rise to slightly larger values of J {Ail.), particularly in the relevant case of cuspy 
DM profiles. 

^This yields a further factor 0.5 in Eq. ((2j for the case of HESS. 
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In this Section we take the opportunity to analyze the recently reported 2004 HESS data in terms of 
DM annihilations, adopting a fully model-independent strategy, which we also apply, for comparison, 
to the HESS 2003 and to the Cangaroo-H results. The purpose of this study is two-fold: first, we 
shall assess the preferred range for the DM masses and pair annihilation cross sections in the top- 
down interpretation of the ACTs data, in the interest of, and as a motivation to, the remainder of 
this paper; second, we wish to extend and generalize the above mentioned analyses of ACTs data, 
privileging a fully numerical and model-independent approach. In this respect, instead of analytically 
approximating the gamma rays spectrum generated by DM pair annihilations, we use the Monte-Carlo 
simulations results from the Pythia code [35] for the gamma rays spectral functions, as implemented 
in the DarkSUSY package [36]. Further, we numerically determine the "&esi spectral functions" , i.e. 
the set of branching ratios minimizing the of a given ACT data set fit, hence providing the most 
conservative confidence level regions for the DM particle mass and annihilation cross section. 

As a first step, we want to pinpoint the pair annihilating DM mass range and final state branching 
ratios pattern favored by the HESS data. Given a set of branching ratios {o/}, such that aj = 1, 
and a WIMP mass m^, one can analytically determine the value of the quantity ((tJ), defined in 
Eq. (jSJ, providing the lowest in the fit to a certain experimental data set {i} with data points 
(4,$^±A$^)), to be 

^ (S^^^F) - [Sr ^ ^^y'- j 

(5) 

As a side remark, in the computation of the x^i ^ well as in Eq. ((21), we take into account the 
finite energy resolution of the HESS detector, which is ~ 15% [20], and consider, for every data point, 
the average WIMP induced gamma rays fiux in each resulting finite energy interval. Eq. Q allows 
to determine the minimal for ^ given WIMP model (with the branching ratios set {o/}), at each 
annihilating particle mass. In order to find the absolute minimum, one needs to find the ^^besf {af} 
set. To this extent, we applied a Monte Carlo technique: we started from a random configuration {a/}o, 
computing, through equation the associated minimal Xq. We then generated a new configuration 
{a/}i, varying one of the coefficients, and re-normalizing the whole set so that X]/^/ ~ ^ (the size 
of the variation was chosen to optimize the convergence to the minimum). The new set {a/}i was 
then accepted with probability 1 if the resulting x? < Xo' with probability exp(— (xi — Xo)/^) 
if Xi > Xo (in order to avoid local minima in the parameter space). The whole procedure was then 
re-iterated, keeping track of the absolute lowest x^ value reached along the resulting "Markov chain" . 
The associated branching ratio set {a/jmin was eventually selected as providing the "6esi spectral 
function" at that particular WIMP mass. 

Besides this fully model-independent approach, we also picked a few representative spectral func- 
tions. In particular, we took the case of a DM particle annihilating with 100% branching ratio into a 
given final state. Based on the spectral features of the gamma ray spectrum (see e.g. Ref. [13,34]), 
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Figure 1: The minimal resulting from fits to the 2003 HESS data (left panel), and to the full 2003+2004 
HESS data set (right panel) on the gamma-rays flux from the Galactic Center, as a function of the annihilating 
Dark Matter particle mass, for different final state channels (see the text for details), including, as a thick black 
line, the "best spectral function" case. The two horizontal lines in each panel indicate the 90% and 95% C.L. 
exclusion limits corresponding to the set of data under consideration. 

we singled out three significant final states, namely (ordering from the softer to the harder spectrum) 
a heavy quark-antiquark pair (for definiteness, bh), a pair of massive gauge bosons (Vl^"*"!^"), and a 
pair of taus {t~^t~), including in the latter case the possibility of a photon in the final state. We 
also addressed the case of Kaluza-Klein (KK) DM, in the context of Universal Extra Dimensions. We 
included in our analysis the internal bremsstrahlung effect discussed in Ref . [30] , and we made use the 
branching ratios obtained in Ref. [37]. The resulting spectrum is quite hard, and significantly different 
from those mentioned above, mainly due to the helicity unsuppressed leptonic final states contributing 
to the process 

5(1)^(1) ^ /+/-^. 

We show our results, for the three final states hb, W~^W~, t~^t~, for the case of KK DM and 
for the "6est spectral function" , in Fig. ^ for the the 2003 HESS data (left panel) and for the full 
HESS data set (2003+2004 data, right panel). We notice that WIMPs purely annihilating in bb or 
in W~^W~ , and KK DM particles give, in different ranges of masses, statistically acceptable fits to 
the 2003 HESS data: this means, in particular, that the spectral features of the annihilating particle 
were not tightly constrained by the 2003 HESS data. Including the new 2004 HESS data, instead, 
statistically rules out, with a high confidence level, both the cases of a WIMP purely annihilating into 
gauge bosons (a final state which applies to many supersymmetric models, including the case of pure 
higgsinos and winos) or quark pairs, and the KK DM scenario [30]. 
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In order to improve the fit, we resorted to a mixed final state composed by a branching ratio x into 
T~^T~ and 1 — X into bb. The resulting spectral function features a harder spectra at E^/m^ 1, and 
improves the fit at the larger energies probed by the 2004 HESS campaign. Further, this particular final 
state is well motivated, in the framework of super symmetric DM. At large tan /3, bino-like neutralinos 
tend to have a final state pattern as that considered above, with a naively estimated relative weight, 
on the basis of color factors only, and neglecting the fermions and sfermions mass effects, x ~ 0.25 = 
1/(3 + 1). In this case, a region allowed at 90% C.L. appears for the mixed final state case. 

We also considered the internal bremsstrahlung of W pairs discussed in Ref. [38], which gives 
a further hard component to the spectral function of the W~^W~ final state at E^/m^ 1. The 
resulting minimal for a WIMP annihilating with branching ratio 1 into charged gauge bosons {e.g. 
a wino-like neutralino) is typically lower than what shown in Fig. ^ but it is a rapidly oscillating 
function of the mass. This simply depends on where the harder end of the spectral function (at 
E^/ < m^) is located with respect to the HESS data points. In any case, we find that the minimal x^'s 
per degree of freedom are always larger than ~ 50/30, ruling out this scenario at 99% C.L.. 

Resorting, finally, to the "&est spectral functions", the relevant, most conservative mass range 
allowed at 90% C.L. reads 

6 < m;^/TeV < 30. (6) 

The WIMP models providing the best fit to the HESS data set are found to feature spectral functions 
given by suitable mixtures of soft and hard channels. On average, we find that the total branching 
ratio into quark pairs ("soft channels") lies between 40% 70%, that into gauge bosons is always 
less than 20% and that into "hard" channels (e.g. charged lepton pairs) lies between 20% 40%, the 
fluctuations in these figures being rather significant, depending on the WIMP mass. 

To have a more exhaustive picture, we show in Fig. I^lthe most conservative (i.e. computed in 
the best spectral function approach) C.L. regions for the Cangaroo-II and HESS data (we separately 
analyze here the 2003 and 2004 data) , on the full plane spanned by the annihilating WIMP mass versus 
the quantity (crJ) defined in Eq. Q. We adopt here, at every m^, the same procedure of Monte Carlo 

minimization described above, dropping, naturally, the use of Eq. ©. As discussed above, the 
2004 HESS data significantly shrank the allowed region, which, quite remarkably, is fully compatible 
with the best fit region of the 2003 data. The different spectral indexes found by the Cangaroo-II and 
HESS collaborations explain the unsimilar shape of the favored (m^, {a J)) regions (particularly in the 
way the best (crJ) values vary with m^), which are clearly incompatible among each other. 

The displayed confidence level contours determine the ranges of (cJ) allowed at a given C.L., and 
hence enable to nail down the relevant intervals on the physical (J(Ail), {av)) plane. This is done 
in Fig. 13 where we jointly analyze the full 2003+2004 HESS data set (the analogous of Fig. [21 is 
shown in the small panel in the upper right). The red arrows on the x-axis indicate the value of the 
quantity J (AO,), defined in Eq. (jlj, computed for the relevant solid angle AJIhessj foi^ four different 
halo model profiles, ranging from the cored Burkert profile [39] to the extremely cuspy adiabatic 
contraction [40] of the N03 profile [41] (ACN03 below) (for the definition of the halo models under 
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Figure 2: Iso-confidence-level contours of "best spectral functions" fits to the Cangaroo-II and to the 2003 and 
2004 HESS data, in the plane defined by the annihilating particle mass and by the quantity (cfJ), defined in 
Eq. fi. 

consideration here see also Ref. [42] and Ref. [36]). The upper x-axis indicates the so-called "5oost 
factor'\ defined as the deviation of J (AO,) from its value for the Navarro Frenk and White (NFW) 
profile, i.e. J( Afi)/ Jnfw^^ [30]. At 90% C.L., and focussing on the two cuspiest profiles under 
consideration, respectively the NFW profile and the ACN03 profile, we find that the relevant intervals 
for the DM annihilation cross section, in units of cm^s~^ and considering the full 2003+2004 HESS 
data set, read: 

2.5 X 10"^^ < {av) < 1.1 X 10"^^ (NFW) 3.0 x 10~^^ < (av) < 1.3 x lO'^^ (ACN03) (7) 

Assuming the naive relation between a WIMP pair annihilation cross section and its relic abundance 

(av) w 3 X lO"^'^ cm^s"V^^x^^ (8) 

one would conclude, from the above listed ranges, that a WIMP providing the WMAP relic abundance 
indicated in Eq. Q could yield a large enough gamma ray flux to explain the ACT data only if an 
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Figure 3: Confidence level contours in the {{crv) , J) plane, for the full HESS data set. The upper x-axis indicates 
the "boost factor" , with respect to the NFW halo profile, i.e. J/Jnfw- The red arrows indicate the values of J 
for a sample of viable Milky Ways halo profiles. In the upper right panel we show the confidence level regions 
on the {m^, i'^J)) plane for the 2003+2004 HESS data, in the "best spectral functions" approach. 

extremely cuspy profile, or a boost factor of the order of 10^ is hypothesized for the DM density in 
the galactic center^. We will show below that the above mentioned relation © can badly fail, and 
that not only can very massive (0(10 -+ 100) TeV) neutralinos produce a relic abundance in the range 
of Eq. but also that, even at those large masses, supersymmetric models allow for (av) values 
implying significantly less cuspy DM profiles than what the estimate above would suggest. 

3 Gravity and anomaly mediated SUSY breaking models 

The phenomenological investigation of supersymmetric extensions of the Standard Model faces the 

ineluctable necessity of handling a huge parameter space, a consequence of our ignorance of the mech- 

^If < ^^CDM, and the dark matter is made up of other components besides the neutralino, the neutralino number 
density squared, relevant for the photon fluxes, would be naturaUy rescaled by a factor (Qchi/^lcDM)'^ , and the needed 
cross sections or boost factors to explain the ACT data would be even larger. 
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anism of supersymmetry breaking (see Ref. [15] for a recent review of the soft supersymmetry breaking 
lagrangian). A way out of this conundrum is to assume that a particular physical process dominates 
the coupling of the hidden sector, where supersymmetry is broken, to the visible sector. A (partial) 
list of such processes includes gravity, (super- Weyl) anomaly and gauge mediated supersymmetry 
breaking models [15,43]. A set of further assumptions, motivated on theoretical or phenomenological 
grounds, has lead to the definition of a few minimal frameworks, featuring a manageable parameter 
space (see e.g. Ref. [44]). The question of how well these minimal setups faithfully reproduce the 
phenomenological features of the general MSSM [45], or of the subset of MSSM models sharing the 
same supersymmetry breaking mechanism, is both long-standing and, in many respects, not totally 
understood. However, in view of the great wealth of phenomenological studies focusing on minimal 
models, it is certainly worthwhile to consider them as benchmark setups, and to start our present 
investigation from them. 

A very popular framework which assumes a gravity mediated supersymmetry breaking scenario is 
that of minimal supergravity (hereafter mSUGRA, see Ref. [46]). Universal scalar and gaugino soft 
breaking masses (respectively niQ and M1/2) plus a common trilinear scalar coupling (^o); all of them 
defined at the grand unification (GUT) scale, and the requirement of successful radiative electro-weak 
symmetry breaking (REWSB) reduce the parameter space of mSUGRA to four continuous parameters 
plus one sign, viz., 

mo, Ml/2, ^0, tan/3, sign(/i) (9) 

where tan f5 indicates the ratio of the Higgs field vacuum expectation values and sign(;u) stands for 
the sign of the supersymmetric fj, term. 

The minimal anomaly mediated supersymmetry breaking (mAMSB) scenario is motivated by the 
possible dominance of supersymmetry breaking contributions originating in the super- Weyl anomaly, 
which are always present when supersymmetry is broken [47]. When the supersymmetry breaking 
and the visible sector reside on different branes, sufficiently separated in a higher dimensional space, 
gravity contributions can in fact be strongly suppressed [48]. In this case, the resulting soft parameters 
are UV insensitive, and can be expressed in terms of low energy entries, such as the Yukawa and gauge 
couplings and the gravitino mass, Foi' instance, the gaugino spectrum is given by 

M, = ^m3/2, (10) 
9i 

where Pg. indicate the beta functions of the gi coupling, i = 1,2, 3. This yields a specific relationship 
in the low energy gaugino mass spectrum, namely 

Ml : M2 : M3 = 2.8 : 1 : 7.1, (11) 

which reverse the mSUGRA hierarchy between the bino and wino mass terms {M2/M1 ~ 2). A 
further universal scalar mass mg is postulated to cure tachyonic sfermions, and radiative electroweak 
symmetry breaking is required, thus yielding an overall parameter space consisting of the following 
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set [49,50]: 



tan/3, 7713/2, mo, sign(^). 



(12) 



3.1 Maximal neutralino masses in Minimal models: mSUGRA &: mAMSB 

Can minimal, benchmark supersymmetric models provide adequate particle physics setups to explain 
the ACT data as analyzed in Sec. As a first step, we need to assess whether the maximal neu- 
tralino masses compatible with the observed DM abundance are large enough to fall within the ranges 
indicated in Fig. ^ Secondly, those large neutralino mass models must feature a neutralino pair 
annihilation cross section compatible with what we nailed down in Sec. [21 

The composition of the lightest neutralino in mSUGRA is dictated by the renormalization group 
equations (RGE) of the gaugino masses from the GUT down to the electro- weak scale, giving, at that 
latter scale, the approximate ratios M2/M1 ~ 2 and M3/M1 ~ 6, and by the absolute value of /i, as 
determined by the REWSB conditions. Since in most of the mSUGRA parameter space |/^| Mi, the 
lightest neutralino is often found to be an almost purely bino-like state. Unless specific mechanisms 
(which will be shortly reviewed below) are in place, pure binos tend to have a thermal relic abundance 
well above the range of Eq. (Q). As a consequence, the parameter space of mSUGRA compatible with 
a viable DM candidate is limited to a handful of special regions, which we list below: 

1. A hulk region at low values of mo and Mij2^ where the t-channel sfermions exchange channel is 
enhanced by relatively light sfermions [51]; 

2. A stau coannihilation region, in the region at low mo where m.^^ ^ m^ [52]; 

3. A resonant annihilation channel through s-channel heavy Higgses {funnel region ), occurring at 
large tan/3, when m^ ~ myi/2, A being the CP-odd MSSM Higgs [53]; 

4. The hyperbolic branch/focus point (HB/FP) region, where large values of mo drive the ^ term to 
low values through REWSB, enhancing the higgsino fraction of the lightest neutralino [54,55]; 

5. A stop coannihilation region, at large values of Aq and low values of mo and M1/2 [56-58]; 

6. A light Higgs resonant annihilation corridor, occurring when ~ nnh/2 [59]. 

It is not difficult to understand that the largest neutralino masses cannot occur in regions 1., 5. 

and 6., where the universal gaugino mass M1/2 is bound to assume low values, well below the TeV 

(and hence the neutralino mass, which cannot be larger, in those regions, than ~ 0.4 x Mi/2)- On the 

other hand, the maximal neutralino mass in the remaining regions can be determined by (a.) setting 

the relic abundance to its maximal value in Eq. and (6.) maximizing the efficiency of the relic 

density suppression mechanism in place. In particular, this amounts to pick points at m^^ = in 

region 2., points at = m^/2 in region 3.^, and pure higgsinos in region 4. [54,55]. 

■^Since the freeze-out process occurs at finite temperature, the minimal rehc abundance in the funnel region strictly 
occurs at 7^ mA/'2\ however, since the induced variation in the maximal mass is small, and since we want to maximize 

also the neutralino pair annihilation cross section at T = 0, we take here the simpler condition = myi/2. 
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Figure 4: (Left panel): Points, in the mSUGRA parameter space at Aq ^ and at = m^-^, giving Q,^h^ = 
0.13, i.e. the largest possible neutralino masses in the stau coannihilation region, for positive (solid blue line) 
and negative (dashed red line) fi. (Right panel): Points, in the mSUGRA parameter space at Aa ~ and at 
= mAl'^, giving ft^h^ = 0.13, i.e. the largest possible neutralino masses in the funnel region, for positive 
(solid blue line) and negative (dashed red line) /i. 

We show in Fig. ^ the resulting maximal neutralino masses in region 2. (left panel) and region 3. 
(right panel), for positive (solid blue line) and negative (red dashed lines) values of sign(;u). The dashed 
line in the left panel terminates where REWSB is no longer possible, while the lines in the right panel 
start where = m^^ . The peaks in the left panel correspond to the overlap of stau coannihilations 
and of resonant annihilations, and give the largest neutralino masses in these two regions, < 900 
GeV. In this analysis, we set Aq = 0; a non-zero value for the trilinear scalar coupling affects the 
left-right mixing in the lightest stau, and hence the neutralino relic abundance. We checked that this 
effect amounts to a factor well within ±10%, with larger relic abundances at large and positive Aq 
and smaller relic abundances at large and negative Aq. Since the relic abundance of the degenerate 
neutralino-stau system scales as m^, taking into account the effect of Aq in Fig. |2 amounts to 
allow an uncertainty factor of ±0(5%). 

As a side comment, the points in Fig. [IJ left, where ruy. ~ m:^^, have been shown, in Ref. [60], to 
play a special role in the suppression of unwanted small scale (sub-galactic) structures, in virtue of the 
effects of a long-lived stau in structure formation. Regarding region 4., the reader is referred to the 
model-independent analysis of Ref. [55]; the maximal neutralino mass in the focus point region ranges 
between 1150 and 1170 GeV, slightly varying with tan/?. As a by-product, we also conclude that 
within mSUGRA, the largest neutralino mass compatible with a thermal relic abundance consistent 
with the WMAP result is less than 1.2 TeV. 
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Figure 5: The pair annihilation cross sections of the largest mass neutralinos in the mSUGRA and mAMSB 
models giving a WMAP neutralino relic density. The shaded areas correspond to the 90% C.L. contours for 
the Cangaroo-II data (left, gauge bosons final state case) and for the HESS data (right, best spectral function 
case, the gauge boson final state being ruled out, see Fig.^ right), for a NFW profile (upper shaded areas) and 
for the adiabatically contracted N03 profile (lower shaded areas). The upper (lower) dotted lines correspond to 
extrapolations of the cross section of mAMSB (HB/FP region of mSUGRA) neutralinos giving a relic abundance 
outside the 2-cr WMAP range. 

In the case of the mAMSB model, the parameter space is much more homogeneous [48-50], and the 
Hghtest neutrahno is almost everywhere wino-like to a high degree of "purity" . Since winos efficiently 
annihilate into gauge bosons, and undergo coannihilations with the lightest chargino, with a small 
spread in the relative mass splitting, provided the lightest neutralino is the LSP its relic abundance 
smoothly depends on its mass, according to the functional relation [61] 

Vt^h^ - c- ( ^^^y )^. 0.0225 < c < 0.0255, 7 -1.90 ^1.92 (13) 

the spread in the coefficients being motivated by the details of the model particles spectrum. As a 
consequence, it is easy to figure out the range of neutralino masses falling within the WMAP range. 
A numerical scan of the mAMSB parameter space indicates that the relic abundance of mAMSB 
neutralinos falls within the range of Eq. for neutralino masses 2.0 TeV < < 2.6 TeV. 

We can already draw, at this point, the conclusion that minimal models do not give large enough 
neutralino masses to explain the HESS data (see ©). In order to understand whether these models 
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feature a large enough pair annihilation cross section to explain, at least, the Cangaroo-II data, we 
map the mSUGRA regions corresponding to the maximal neutralino masses, and the mAMSB points 
with WMAP relic abundances, on the {m^, {cv)) plane. The neutralino pair annihilation cross sections 
have been computed, here and in the remainder of this paper, with the DarkSUSY package [36]. We 
shade in light blue the 90% C.L. contours corresponding to fits to the Cangaroo-II data with a 100% 
branching ratio into gauge bosons pairs, the case relevant for the largest neutralino masses under 
consideration here. We also include, for completeness, the best spectral function 90% C.L. region 
corresponding to the (2003+2004) HESS data. The two upper contours refer to the NFW profile, 
while the lower ones to the ACN03 profile. We also extrapolated the mAMSB (pure wino) and FP 
(pure higgsino) lines to models featuring a lower (left portions of the dotted lines) and higher (right) 
thermal relic abundance. We recall that low thermal relic density models can provide all the required 
DM in the context of non-standard cosmological scenarios (including a quintessential energy density 
term [62], Brans-Dicke- Jordan theories [63,64] and anisotropic cosmologies [63,65]) or of non-thermal 
production [66]. We henceforth conclude that 

• the mSUGRA model is unfit to explain ACT data, and that 

• the mAMSB model can explain the Cangaroo-II data only, and with (1) an extremely cuspy DM 
profile as well as (2) some kind of cosmological relic density enhancement mechanism. 

3.2 Minimal models with non-universal Higgs masses 

A common and reasonable approach to fill the gap between minimal models and the general MSSM 
is to relax some of the assumptions of universality, possibly following GUT-motivated guidelines. 
In the present framework, relaxing universality in the sfermions sector [67] would not greatly help, 
since while larger neutralino masses could be made compatible with the DM abundance, the pair 
annihilation cross section would be hardly affected. The two remaining options concern the gaugino 
and the Higgs sector of soft supersymmetry breaking masses. Both possibilities are well motivated 
by supersymmetric GUTs, and have been widely investigated in the literature [68,69]. We pick here 
the option of non-universality in the Higgs sector, which has been recently analyzed in great detail in 
Ref. [70] for the case of mSUGRA. We consider here, for the first time, general soft supersymmetry 
breaking Higgs mass extensions to the mAMSB as well. It will be shown, in the following Sec. 21 
that resorting to what we will hereafter dub as non-universal Higgs masses (NUHM) models actually 
allows to access the largest possible neutralino pair annihilation cross sections of the general MSSM. 

Considering general soft breaking mass terms nin^ and mn^ for the Higgs doublet superfields Hu 
and Hfi amounts, through the conditions of REWSB, to extend the minimal model parameter spaces of 
Eq. @ and H12|) to the extra weak-scale parameters ^ and rriA [69]. Stretching the neutralino masses 
to the largest possible values compatible with a WMAP thermal relic abundance squeezes the NUHM 
parameter space to the regions where resonant neutralino annihilation occurs (m^ ~ m^/2). Further, 
since the neutralino couplings to the Higgses scale as the product of the gaugino and of the higgsino 
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Figure 6: (Left): Curves of iso-relic abundance for neutralinos in the mSUGRA model with non-universal Higgs 
masses, in the {niA, fJ') plane. The green shaded area corresponds to the 2-a WMAP range. (Right): Curves of 
iso-cross section for neutralinos in the mSUGRA model with non-universal Higgs masses, in the {mA,fi) plane. 
The green shaded area corresponds to the 2-a WMAP range, while the yellow shaded area indicates the reach of 
stage-3, ton-sized direct detection experiments. 

fractions, the lowest relic abundances are found in the regions where a maximal higgsino-gaugino 
mixing occurs. 

As sample illustrative cases, we considered, in Fig. El and EJ two NUHM parameter space slices 
along the {mA, fJ,) plane, fixing 

mo = 25 TeV, M1/2 = H TeV, Aq = 0, tan/3 = 40, sign(/i) > (14) 

for the NUHM mSUGRA case (Fig.©, and 

tan/? = 40, m3/2 = 4600 TeV, mo = 25 TeV, sign(/i) > 0. (15) 

for the NUHM mAMSB case. The neutralino masses range, in Fig. El from 4.5 to 5.2 TeV, and in Fig.[7| 
between 11.5 and 12.1 TeV, well above the quoted maximal mSUGRA and mAMSB neutralino masses. 
The left parts of the Figures show the contours at constant relic abundance, and feature a shaded green 
area corresponding to the DM abundance range of Eq. . The neutralino mass can be read out of the 
plots as the minimum between Adi ~ 5.2 TeV (resp. M2 — 12.1 TeV) and /i (on the y-axis). In both 
cases, we see that the lowest relic abundance corresponds to the resonant ~ itia/^ condition, and 
to the maximal higgsino-gaugino mixing region, viz. Mi^2 — y^- The right panels indicate the contours 
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Figure 7: As in Fig. but for neutralinos in the mAMSB model with general soft supersymmetry breaking 
Higgs masses. 

of constant (crv), and the projected reach of future, Stage-3 ton-sized direct detection experiments, 
such as Xenon 1-t [71] (yehow shaded regions). In both cases, the maximal pair annihilation cross 
section region and the WMAP relic abundance areas largely fall within the reach of direct detection 
experiments, despite the corresponding neutralino masses being huge, respectively around 5 and 12 
TeV. Though remarkable, this is not surprising, since the scalar neutralino scattering cross section 
off protons is mediated, besides squark exchanges, by t-channel CP-even Higgses exchanges. In the 
region of large higgsino-gaugino mixing these latter channels have enhanced couplings, analogously to 
what happens for the neutralino pair annihilation resonant CP-odd Higgs exchange. This points to 
the conclusion that, in these models, although multi- TeV neutralinos will be (kinematically) beyond the 
reach of future accelerators, a crossing symmetry between neutralino pair annihilation processes and 
neutralino scattering off matter, enforced by the requirement of a sufficiently low relic abundance, will 
imply detectability at future direct detection experiments. Further, comparing the right panels of Fig. |21 
and d with Fig. El we notice that the range of masses and annihilation cross sections nicely falls within 
the preferred HESS range for the NUHM mAMSB model, and (extrapolating {av) ~ const x m"'^) 
within the Cangaroo-II range for the NUHM mSUGRA model. Concerning the spectral functions 
resulting from NUHM mAMSB models at large neutralino masses, we notice that the sufficiently large 
branching ratios into tau pairs, particularly at large tan/3, yield a rather hard spectrum at P^ — > m^, 
nicely compatible with the 2004 HESS data. 
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Figure 8: (Left): The thermally averaged effective cross section ((Toffw)(T), at a neutralino mass of 1 TeV, for 
four models (respectively a mAMSB wino-like neutralino, a mSUGRA neutralino lying in the HB/FP region, 
and two resonantly annihilating models at mi ~ fi (mixed bino-higgsino) and m2 ~ fJ-, such that ~ mA/2.). 
(Right): cross sections of resonantly annihilating neutralinos in the (m^, (cry)) plane. The green strips corre- 
spond to models giving a thermal neutralino relic abundance within the 2-a WMAP range. The light blue and 
red shaded areas refer to the 90% C.L. regions for the Cangaroo-II data and for the 2003+2004 HESS data, 
for a NFW profile (two upper shaded areas) and for the adiabatically contracted N03 profile (two lower shaded 
areas). We assume here that neutralinos pair annihilate with BR(bb)=70% and BR(t^t^ )=30%. 

4 The largest annihilation cross sections 

What is the largest possible neutralino pair annihilation cross section, at a given neutralino mass, in 
the minimal supersymmetric extension of the standard model? 

Despite the MSSM parameter space being so large, one can make general statements regarding 
the supersymmetric configurations featuring the largest possible values of (cv). For instance, given a 
certain particle spectrum, since binos do not couple to gauge bosons, and feature a suppressed s-wave 
amplitude, the neutralino LSP composition of the models giving the largest possible annihilation cross 
sections, banning special cases, as resonant annihilations, will be mostly wino or higgsino-like, or a 
mixture of the two. We saw above that the sfermion spectrum is not critical for the wino/higgsino- 
like pair annihilation cross section: the chargino or neutralino t-channel exchange amplitudes always 
dominate over sfermion exchange ones. The remaining possibility therefore pertains to the Higgs 
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sector, and to the occurrence of s-channel resonances (see Ref. [72] for an accurate account of the 
treatment of resonant annihilation channels). The only MSSM particle giving a non- vanishing s-wave 
amplitude is the C-P-odd A Higgs boson [14]. Needless to say, the maximal resonant effect will be 
obtained at niy. = and, as discussed above, at maximal higgsino-gaugino mixing. In order to 

compare the pair annihilation cross section of resonantly annihilating, maximally gaugino-higgsino 
mixed neutralinos with the standard effective cross section of winos and higgsinos, we show in Fig. |S1 
left panel, the thermally averaged {aef[v)iT) (see Ref. [73] for the definition of aes, i-e. the effective 
cross section including properly weighed coannihilation effects), at a neutralino mass of 1 TeV, for 
four models (respectively a mAMSB wino-like neutralino, a mSUGRA neutralino lying in the HB/FP 
region, and two resonantly annihilating models at mi = fi (mixed bino-higgsino) and m2 = /x, with 
= myi/2). Evidently, models with an open resonant annihilation mode give the largest cross 
sections at T = 0. 

In order to have an analytical insight, let us consider (av) in the limit of zero relative velocity, at 
= mA/2 

2 I 2" 

{<^v)iT = 0) = ^ E^fPf\9ffA\', /?/ = Jl - 5' (16) 

A f V X 

where / stands for the final state fermions, c/ is a color factor equal to 3 in case of quarks and to 1 in 
the case of leptons, g^^^A and 5//A respectively indicate the neutralino pair and fermion-antifermion 
couplings to the CP-odd Higgs boson A, and Ta is the A total width. Disregarding QCD corrections, 
we can write 

and 

9xxA = {92Z12 - giZii) (Zi4 cos/3 - Zi3 sin/?) , (19) 
using the following convention for the lightest neutralino composition: 

X = ZuB + Zi2W^ + ZisH^ + Zi4i?°. (20) 

Neglecting the two lighter fermion generations, and assuming a heavy enough neutralino (namely, 
rriy- ^ rrit), in order to maximize Eq. (|16() . in view of the mild variations with tan/? of Eq H19|) . we 
need to find the maximum, as a function of tan /3, of 

((3mg + m^)tan2/3 + 3m?/tan2/3)"^ (21) 

which is found to occur at 

tan/3 ^mt/rrib 6 ^ 7. (22) 
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Plugging the value of the minimum of Ta into Eq. H16|) . we give the following estimate, 




(23) 



the subscripts 1 and 2 referring to the case mi < m2 and mi > m2 respectively (i.e. to a mixed bino 



Eq. H23j) shows that the maximal annihilation cross section occurs for maximally mixed (last factor) 
wino-higgsino neutralinos (since 52 > <?i)- The maximal bino-higgsino annihilation cross section will 
be suppressed, according to Eq. (j^Hj), by a factor (51/^2)^ = tan^ 0.3. 

The above predictions are nicely confirmed by our numerical results, showed in the right panel of 
Fig.lSlon the (m^, {av)) plane. We picked three ^^pure" and two ^^maximally mixed' MSSM neutralinos, 
setting the (largely unimportant) common sfermion mass to m^, = 10 x m^, tua = 2 x m^, and the 
following relations for the gaugino masses and for the term, at the weak scale, 

fi = 10 X mi, m2 = 2 X mi, = 6 x mi (Bino — like, mSUGRA relations) 

mi = 10 X ^, m2 = 2 X mi, = 6 x mi (Higgsino — like, mSUGRA relations) 

/i = 10 X m2, mi = 3.3 x m2, = —8 x m2 (Wino — like, mAMSB relations) 

/i = mi, m2 = 2 X mi, = 6 x mi (Mixed Bino — Higgsino, mSUGRA relations) 

/i = m2, mi = 3.3 x m2, ma = —8 x m2 (Mixed Wino — Higgsino, mAMSB relations) 

For all models, the value of tan /3 was chosen, at each m^ in order to numerically maximize {av) , and 
found to be always close to our estimate above, tan/3 w 7, except for < 4 x m^, where even lower 
values of tan /? gave larger cross sections. In the Figure, we indicate with a thick green line the range 
of masses giving a thermal relic abundance within the WMAP range of Eq. As a side remark, 
notice how badly the naive relation quoted in Eq. Q fails for resonantly annihilating models. The 
light blue and red shaded areas correspond to the 90% C.L. fit countours for the Cangaroo-II data and 
for the 2003+2004 HESS data (computed for a final state with BR(66)=70% and BR(t+t~)=30%, 
indicative of the final state pattern of resonantly annihilating neutralinos). We notice that the HESS 
data can have a supersymmetric interpretation with a model producing the correct thermal neutralino 
relic abundance, and with a boost factor o/O(10). 

Assuming a GUT unification relation for the gaugino masses, a mixed wino-higgsino state cannot 
be realized. Henceforth we conclude this Section quoting the two relevant {av) theoretical upper limits 
obtained: 



or wino LSP), and Zh = Zf^, Zg^ ^ = Zf^ Inserting the numerical values, we obtain 



{av){T = 0)^ £ii±MGeV-2 ~ lO'^s ^ 10-22 (^^\ ' cmV 
mt V 1 TeV / 



1 
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Figure 9: The (av) contributing to the continuum gamma-rays yield from neutralino pair- annihilations including 
non-perturbative electro-weak effects [16], for a wino-like (solid blue line) and for a higgsino-like (dashed red line) 
lightest neutralino. The chargino-neutralino mass splitting was computed including loop effects from custodial 
SU (2) -breaking in the gauge boson sector, i.e. in the limit of heavy sfermions [75]. We also indicate the 90% 
C.L. regions for the Cangaroo-II data set fit, for neutralinos annihilating into gauge bosons. We recall that the 
pure gauge bosons final state is excluded by the 2004 HESS data (see Sec.\^. 

We stress that the results quoted above apply to neutralino masses uiy. > 100 GeV, since all the 
models we considered feature a quasi-degenerate chargino, whose mass is bounded to be larger than 
approximately 103 GeV [74]. 

4.1 The role of non-perturbative EW effects 

Recently, it has been noticed that SU{2)l non-singlet neutralinos (i.e. wino or higgsino-like states), 
featuring a quasi-degenerate chargino partner, with mass in the TeV range or above, undergo, at 
small velocities, large enhancements in their pair annihilation cross section [16,76,77]. This is due 
to non-perturbative effects, leading to the formation of multiple neutralinos-charginos "bound states" 
which can resonantly contribute to the annihilation cross section in the non-relativistic limit [16]. In 
Ref. [16,76] a non-relativistic effective action for SU{2)l doublets (Higgsinos) and triplets (Winos) 
was derived, leading to a numerical evaluation of the pair annihilation cross section of neutralinos into 
final state gauge bosons W. This non-perturbative contribution was never taken into account before. 
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and a convenient numerical fitting formula for (cr?;)v\/') as a function of the neutralino mass and of its 
mass splitting 5m with its chargino partner, was provided as well [16]. Evidently, non-perturbative 
electroweak effects go beyond the treatment outlined above, and must therefore be analyzed to assess 
if and where the theoretical upper bounds on (av) quoted in Eq. H25|) . H26() are violated. 

In the present context, we sum over the gauge boson final states contributing to the gamma rays 
continuum, namely W^W~ and ZZ. We assume, for definiteness, and to get a model-independent 
result, that the dominant loop corrections entering 5m come dominantly from the gauge bosons loop 
contributions [75, 78], which only depend on the neutralino mass. We show our results in Fig. |^ 
The thick solid blue and the red dashed lines reproduce the numerical fitting formula of Eq. (61) 
hep-ph/ 0412403j /2 within their range of validity. The black line indicates the ("perturbative") 
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upper limit quoted above in Eq. Evidently, when hits those special values (depending on 

the chargino-neutralino splitting) corresponding to "zero energy" resonances, {av) locally exceeds the 
perturbative result of Eq. (|25|) . However, this happens in very limited neutralino mass ranges. Further, 
it would be very difficult to reliably predict the value of {av){m^)^ which is a very steeply varying 
function, close to the resonances. The non-perturbative electro- weak resonance scenario appears in any 
case unfit to explain the HESS data. In fact, as pointed out in Sec.|21 the resulting pure gauge bosons 
final state spectral function for the gamma rays continuum (even including the internal bremsstrahlung 
effects of Ref. [38]) is statistically excluded by the analysis of the 2004 HESS data. 

5 The largest masses: non-perturbative QCD effects and gluino 
coannihilations 

In the previous Section we found that neutralinos with the largest possible pair annihilation cross 
section produce a thermal relic abundance compatible with the WMAP upper bound quoted in Eq. 
for masses my^ < 12.5 TeV. It is well known, however, that when a particle exists that can pair 
annihilate [^^co- annihilated^) with the lightest neutralino and that is close in mass with it, so that 
their freeze-out processes overlap and interfere in the early Universe, the final thermal relic abundance 
of the neutralino can be largely affected [79]. We therefore expect that the largest neutralino mass 
compatible with a thermal production of DM can be found in models where coannihilation processes 
are active. The s-wave unitarity limit on a thermally produced WIMP mass derived in Ref. [80], 
suitably rescaled with the result of Ref. [1] quoted in Eq. leads to a putative upper bound of 
my' < 120 TeV. However, the s-wave unitarity argument does not always apply (a counter example is 
given for instance by resonantly annihilating WIMPs, as those considered in the previous Section). We 
wish to assess here which is the maximal theoretically allowed neutralino mass in the general MSSM, 
and whether or not the unitarity bound quoted above applies. 

As a rule of thumb, if the neutralino pair annihilation cross section is smaller than the coannihilating 
partner pair annihilation cross section times the ratio of the number of internal degrees of freedom 



21 



of the neutralino over that of the coannihilating partner^ the neutrahno reUc abundance in presence 
of coannihilations will be reduced. Again as a rule of thumb, the Maxwell-Boltzmann equilibrium 
distribution fixes the relative "weight" of the coannihilation contribution to the effective cross section 
determining the neutralino relic abundance [79] as scaling oc exp(— Am/Tf.o.), where Am is the mass 
splitting between the neutralino and the coannihilating partner, and Tf o is the neutralino freeze-out 
temperature. 

Many minimal supersymmetry breaking models include the possibility of coannihilation processes, 
for instance the mSUGRA and the mAMSB scenarios discussed above. It goes without saying that 
the largest the pair annihilation cross section of the coannihilating partner, the largest the maximal 
possible coannihilation effects. In Ref. [81] a model- independent analysis showed that, not surprisingly, 
the "strongest" coannihilating partners in the MSSM are strongly interacting particles (SIPs), the 
squarks and the gluino. Stop coannihilations were considered previously in Ref. [56-58]. The analysis 
of Ref. [81] was based on the micrOMEGAs code for the relic density computation [82], which includes 
the perturbative cross section only for SIPs. Already at the perturbative level, gluinos were found to 
be by far more efficiently coannihilating than squarks. It was also noticed, however, that special MSSM 
realizations, including the occurrence of multiple squark coannihilations or of resonantly annihilating 
coannihilation partners, could yield net effects on the neutralino relic abundance even beyond that 
of the (perturbative) gluino coannihilations. To quantitatively assess this possibility, we carried out 
an extensive scan of the general MSSM parameter space, along the lines of Ref. [45] (where the 
reader is directed for further details). We find that the above mentioned special, extremely fine-tuned 
models where peculiar squark coannihilations occur can suppress the neutralino relic abundance to 
the appropriate DM density level for neutralinos as heavy as 20-^25 TeV. 

The inclusion of non-perturbative strong interactions effects in SIPs pair annihilations clearly plays 
a crucial role for an accurate neutralino relic abundance computation, and henceforth to nail down 
the maximal neutralino mass compatible with a thermal relic abundance in the CDM range. In 
Ref. [56] it was claimed that a full higher order QCD calculation in the stop (co-)annihilation cross 
sections can amount to corrections around a factor 2 or so with respect to the perturbative result. 
Strong interactions effects are, on the other hand, expected to be much more relevant in the case of 
gluino pair annihilations. With this in mind, we consider below the case of gluino (co-)annihilations, 
motivated by the available results of the extended and accurate analysis of the non-perturbative QCD 
effects in the gluino pair annihilation cross sections carried out in Ref. [83]. 

Supersymmetric scenarios where gluinos can coannihilate with neutralinos include the widely dis- 
cussed class of models featuring non- universal gaugino masses [68]. A string-inspired model which 
favors a scenario where the gluino can be either the LSP or a (possibly coannihilating) next-to-LSP 
is the O-II string model in the limit where supersymmetry breaking is dominated by the universal 
"size" modulus [84]. A recently proposed framework where gluinos can be coannihilating partners of 
a neutralino LSP is that of Split-Supersymmetry, where the superpartners scalar sector is assumed 
*For a more precise condition, including the co-annihilation cross sections contribution, see Ref. [58] 



22 



to be "split" from the fermionic sector [61,85]. Other specific theoreticaUy motivated models were 
discussed in Ref. [45,83]. 

The publicly available numerical package for the computation of the neutralino relic density 
DarkSUSY [36] does not include gluino coannihilations, while, as mentioned above, the micrOMEGAs 
package [82] does not go beyond the leading perturbative cross sections for gluino (co-)annihilations. 
We therefore developed an independent numerical code to account for gluino coannihilations, including 
various non-perturbative gluino pair-annihilation scenarios, and we interfaced it with the {ae{iv){T) 
computation of the DarkSUSY package to compute the neutralino relic abundance, hence automatically 
taking into account all the effects of other coannihilating partners, of resonances and of thresholds. 
We do not include higher order QCD effects in the gluino-neutralino cross section, since we do not 
expect them to be quantitatively relevant. Further, as already observed in Ref. [45], even consid- 
ering a fully perturbative gg annihilation cross section, the gluino-neutralino coannihilation term is 
always sub-dominant in the {a(,ffv){T) computation (it was found there to give at most a per-cent 
contribution). We stress here, however, that gluino-neutralino conversion processes are essential to 
keep gluinos in thermal equilibrium with neutralinos through scattering off relativistic quarks during 
freeze-out. Since those processes are mediated by squark exchanges, this means that squarks cannot 
be exceedingly heavy, not to effectively decouple the gluino and neutralino freeze-out. 

The dynamics of gluino pair annihilations, and the related issue of the relic abundance of a gluino 
LSP, was widely discussed in the literature [83,86]. The perturbative gluino pair annihilation cross 
section as a function of the center-of-mass energy s reads, in the notation of Ref. [83] 

o-p = (T{gg ^ gg) + ^ a(gg qq) (27) 

where 

crigg - gg) = ^ (log ^(21 - e/?^ - 3(3^) - 33/? + 17P^^ , (3 = ^l-imi/s 

a{gg^qq) = ^{3-(3'){3-p'), P=^l-Am^^/s. (28) 

In Ref. [83] it has been claimed that non-perturbative effects, relevant when ^/s ~ 2m^, are expected 
to range between two extreme scenarios. In the first scenario, one only considers the effects of mul- 
tiple gluon exchanges, which, neglecting the logarithmic enhancements due to soft radiation, can be 
parameterized by the (exponentiated form of the) the Sommerfeld enhancement factor, 

where as{Q) is evaluated at the scale of the typical momentum transfer of the exchanged soft gluons, 
Q ~ Pm-g, and the Cgg and Cqq respectively apply to the gluon pair final state cross section and to 
the quark-antiquark pair. 
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Figure 10: (Left): the thermally averaged gluino pair annihilation cross section as a function ofrUg/T (lower 
X axis) and of T (upper x-axis, respectively) , for a gluino mass m-g = 1 TeV, assuming a fully perturbative up, 
a Sommerf eld- enhanced E x ap, and a fully non-perturbative annihilation cross section. We also include the 
s-wave unitarity limit [80] (solid black line) and the cross section of a 1 TeV purely wino-like neutralino (brown 
dot-dashed line). (Right): The effective thermally averaged neutralino pair annihilation cross section (av) (i.e. 
including co-annihilations) as a function ofm^/T (lower x axis) and ofT (upper x-axis), for a mSUGRA bino- 
like neutralino with a 1 TeV mass (tan (3 = 30, toq = ^ > 0, = 0^, without gluino coannihilations 

(solid black line) and including gluino coannihilations. In the latter case we show the lines corresponding to the 
four gluino annihilation cross sections scenarios of the left panel, for gluino-neutralino mass splittings such that 
the resulting final neutralino relic abundance corresponds to the central non-baryonic Dark Matter abundance 
deduced from the WMAP data [1], fl^h'^ ^ 0.113 . 

In the opposite, extreme non-perturbative scenario, gluinos undergo, at small /3's, a transition into 
color singlet bound states. The non-perturbative pair annihilation of those states into vr's is usually 
assumed to be 

aNP = A(3-'/ml, A = 0{1). (30) 

The location of the transition is highly uncertain; in the most extreme scenario considered in Ref. [83], 
the transition was taken to occur when the total gg kinetic energy in the center of mass frame fell 
below a given limit L, with L ~ 0.2 1 GeV. Below the gluino-bound state transition, the scattering 
cross section was assumed to follow the perturbative result in Eq. (|27jl . with either an abrupt or 
a smooth transition to the non-perturbative annihilation cross section (|30|) . The sudden transition 
option was found, in Ref. [83], to give a smaller final gluino relic abundance, and hence we expect, 
with this choice, the maximal possible gluino coannihilation effects. 
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We show the ghiino pah' annihilation cross sections we consider here in Fig. IIUI left panel. There, 
we plot {av) as a function of m/T and of T (lower and upper x-axis), for a gluino mass of 1 TeV. 
We picked two options for the non-perturbative scenario, featuring a transition at a kinetic energy of 
(1) 1 GeV and (2) equal to the vr mass. We also plot the s-wave unitarity limit and the annihilation 
cross section of a 1 TeV mAMSB wino, for comparison. The extreme non-perturbative scenario above 
violates the s-wave unitarity limit on the pair annihilation cross section of a WIMP, as already pointed 
out in Ref. [83], as a consequence of coherent contributions from multiple partial waves. It can be 
noted that at large temperatures all scenarios converge to the perturbative cross section, but at low 
temperatures the differences in the various pair annihilation options can be substantial. 

The importance of including gluino coannihilations for the relic abundance of neutralinos can be 
appreciated in the right panel of Fig. IIUI We consider, there, a 1 TeV bino-like neutralino taken from 
the mSUGRA parameter space, at mo = 5Mi/2) A* > 0, tan/3 = 30 and Aq = 0, and we indicate its 
thermally averaged cross section with a thick black solid line. The relic abundance of that neutralino 
without gluino coannihilations would be Vl^h? ~ 27. We then add a coannihilating gluino, featuring a 
pair annihilation cross section following one of the four schemes of Fig. 1101 left, with a mass such that 
the resulting neutralino relic abundance gives the central WMAP value, Q.^h? ~ 0.113 (this amounts 
to considering gluino-neutralino mass splittings ranging from 7% in the fully perturbative case to 38% 
in the extreme non-perturbative case with transition at a kinetic energy of 1 GeV). We then plot the 
resulting effective neutralino annihilation cross section , as defined in Ref. [73], i.e. the cross section 
including all the properly weighed (co-)annihilation contributions. We see that around the freeze-out 
temperature, the effective annihilation cross section including coannihilations exceeds the neutralino 
pair annihilation cross section by even four orders of magnitudes; also notice that in order not to 
exceedingly suppress the final neutralino relic abundance, in the extreme non-perturbative scenarios 
the mass splitting is much larger, and hence the temperature at which gluino coannihilations effectively 
disappear is also larger. 

To assess the maximal neutralino mass compatible with the WMAP result on the DM abundance 
in the neutralino-gluino coannihilation scenario, we show contour levels of iso-relic abundance in the 
{ruy-, {mg—m^)/m^) plane in Fig.^fl for the four gluino pair annihilation scenarios of Fig.^1 The gray 
shaded area in the left part of each panel indicates gluino masses excluded by Tevatron searches [74] , 
while the green band locates the 95% C.L. range of the CDM abundance as determined by WMAP [1]. 
While in the most conservative non-perturbative scenario (Sommerfeld enhanced cross section, top 
right panel) the maximal neutralino mass is even below 10 TeV, the two extreme non-perturbative 
frameworks allow for neutralino masses well above 100 TeV, without invoking a terribly small mass 
splitting between the gluino and the neutralino mass (the mass splittings in the stau coannihilation 
region are always well below what we consider here). 

In the extremely non-perturbative gluino coannihilation scenarios considered in the two panels at 
the bottom, the violation of the s-wave unitarity limit in gluino pair annihilations (Fig. 1101 left) is 
effectively transferred^ to the neutralino sector. As a consequence, in those scenarios, neutralinos 
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Figure 11: Iso-level curves corresponding to given neutralino relic abundance values, in the plane defined by the 
neutralino mass (x-axis) and by the relative gluino-neutralino mass splitting (y-axis), for four different gluino 
pair annihilation scenarios. The green shaded area corresponds to the 95% C.L. range for fl^h^ , while the area 
shaded in grey, in the right part of each panel, indicates the region in the plane excluded by Tevatron searches 
for gluinos [74]- 

with masses at, or even above 100 TeV produce a sufficiently low thermal relic abundance, provided 
the gluino has a mass within, say, 20% of the neutralino mass. 

As a bottom line, in view of our discussion above, we conclude that taking into account non pertur- 
bative effects in gluino pair annihilations, and depending on the assumed non-perturbative scenario, the 
largest neutralino mass in the MSSM compatible with a WMAP relic abundance ranges from around 
20 TeV to well above 100 TeV. 
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6 Outlook: the largest "supersymmetric factors" for neutralino in- 
direct detection 

We provided in Sec. ^ an analytical and numerical estimate of the theoretical upper bound on the 
neutralino pair annihilation cross section (cu), as a function of its mass. Neutralinos featuring that 
[av) produce a relic abundance in the WMAP range only at very large masses, around 12 TeV. At 
smaller masses, those models do not thermally produce enough relic neutralinos to explain the observed 
DM abundance. In that case, one either assumes that neutralinos only make up for a fraction of DM, 
the rest being composed by some other particles, or relaxes the assumptions that lead to the usual 
thermal relic abundance result. In this latter approach, as already mentioned, a few loopholes have 
been considered so far: neutralinos can be produced non-thermally in the decay of moduli, gravitinos, 
Q-balls, cosmic strings etc. [66], or a modified behavior for the Hubble expansion factor H (driving the 
neutralino freeze-out) might have occurred, affecting the final relic density [63,65]. Various examples of 
modified cosmologies have been proposed so far (including quintessential scenarios [62], Brans-Dicke- 
Jordan cosmologies [63,64], anisotropic cosmologies [63,65] etc.). 

In the most conservative approach, without assuming non-thermal production or cosmological relic 
density enhancement processes, any signal from neutralino pair annihilations will be proportional to 
the quantity 

e{av)/ml, e = min[l,(0x/i2)/(f^CDM/i')min], (ncDM/i')min = 0.095 (31) 

the factor ^ taking into account the fraction of DM composed by neutralinos. In this case, it is highly 
non-trivial to theoretically predict which is the maximal supersymmetric factor {av) / . We thus 
resorted to a very large numerical scan of the MSSM, including all phenomenological constraints 
{e.g. from supersymmetric contributions to BE(6 57), BR(i?s — > precision electro-weak 

observables, supersymmetric particles searches, Higgs searches, etc.). The scan was again performed 
along the lines discussed in Ref. [45], where the reader is directed for further details. As a by-product, 
we also explicitly verified that no MSSM models in our scan violate the theoretical upper on {av) 
quoted in Eq. (|25j) . We show our results in Fig. ^1 for the pair annihilation cross section (left) and for 
the "supersymmetric factor" (right). We also indicate the region excluded by the neutralino- induced 
synthesis of ^Li in the early Universe, as computed in Ref. [87]. The ^Li constraint strictly applies 
only to the dark blue region, as long as the above mentioned non-thermal production or relic density 
enhancement processes are not assumed. 

Most of the WMAP models featuring a large {av) have, again, a resonant annihilation channel, 
with a slightly off-resonance kinematical condition on ruA ~ 2 x m^. While the lightest neutralino 
mass WMAP models are mostly bino-like, a non-trivial higgsino component starts at ruy. ~ 500 GeV, 
giving a significant increase in the largest possible {av) values. The largest pair annihilation cross 
sections and supersymmetric factors we find in the range we explored {m^ > 100 GeV) are listed 
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Figure 12: The regions, in the {m^, (cu)) plane, populated by super symmetric models with non-rescaled cross 
section and arbitrarily low relic abundances (light-blue hatched region) and by models with a WMAP relic abun- 
dance or with a rescaled cross section (i.e. if f2^/i^ < ft^^h'^ then (av) —y (av) x (fl^/n^^)^ ), indicated by 
dark blue squares. The latter were obtained through a random scan over the general, flavor diagonal and CP- 
conserving, MSSM parameter space. The region lying above the red line produces an excess of ^ Li in the early 
Universe, if all DM is composed by neutralinos. (Right): same as for the left panel, but for the "super symmetric 
factor" {av) /m^. 

below, for the convenience of the reader, 

(^^(az;))^^^ ~ 8 X 10"24 cmV^ at ~ 550 GeV (32) 
{i'^{av)/m^^f''''^ ~ 6 X 10~29 cm^s-^GeV-^ at ~ 150 GeV (33) 

A similar analysis, extending to even smaller neutralino masses than those we consider here, was 
recently carried out in Ref. [88]. Although we find, as in [88], that the maxima quoted in 1)32(1 and 
H33() correspond to models with ^ = 1, the absolute maximal value for {^"^ {av))^^'^ we find here 
is more than two orders of magnitude larger than what quoted in [88], even taking into account 
an uncertainty factor of order unity. The reason for this discrepancy can be traced back to either 
the slightly different MSSM parameters over which the scans have been performed, or to numerical 
differences in the evaluation of the pair annihilation cross section, particularly in the proximity of 
resonances, where we find the largest cross sections reported here. On the other hand, Ref. [88] finds 
larger values for the supersymmetric factor than what we quote in ()33|1 . outside the range of masses 
considered here, at <C 100 GeV. As a side remark, we would like to point out that in presence of 
coannihilation processes, and particularly in view of our results presented in Secl^l no direct relations 
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Figure 13: The regions, in the {m^,{av)) plane of super symmetric models with non-rescaled cross section 
(yellow shaded area) and of models with a WMAP relic abundance or with a rescaled cross section (orange 
shaded area) excluded by the Egret data [32] and/or by the 2004 HESS data [21] on the gamma ray flux from 
the center of the Galaxy, for three different halo models. The thin dashed lines indicate the shape of the Egret 
and HESS bounds, for the NEW profile, where one constraint is weaker than the other. We also indicate, again 
with a thin dashed line and for the NEW profile, the bound resulting from the 2003 HESS data. (Right): same 
as for the left panel, but for the "super symmetric factor" (av) /tti^ . 

link the neutralino pair annihilation cross section with its thermal relic abundance. In this respect, it 
is very hard to draw a minimal, guaranteed value for either ^'^{crv) or ^'^{av) /m'^, which in presence 
of a "strong" coannihilating partner can essentially be arbitrarily low. 

The results presented in this Section might be of relevance in the assessment of the scientific 
purposes of experimental devices sensitive to DM annihilation products, or in establishing whether 
the results of a given experiment provide, or not, any constraint on viable supersymmetric models. 
Further, our results can be of some use in order to construct best case scenarios for DM indirect 
detection signals, e.g. in evaluating the theoretically maximal supersymmetric DM induced gamma 
ray or neutrino fluxes from a given astrophysical system. 

As an illustrative worked-out case study, we consider the halo-model dependent bounds put by the 
Egret [32] and by the HESS [18,20] data on the gamma ray flux from the galactic center on the planes in 
Fig. 1121 Given a DM density profile, one can conservatively require that the neutralino annihilations- 
induced gamma ray flux from the Galactic center alone does not exceed the fluxes measured by 
the Egret and by the HESS experiment, in any of the energy bins of the respective data sets. We 
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performed this exercise in Fig. 1131 where we indicate with a solid black line the contours above which 
the parameter space is excluded at 95% C.L. for three halo models: the Burkert profile [39], the NFW 
profile [31] and the ACN03 profile [41,42]. We also indicate, with a thinner dashed line, the shape 
of the regions exclued by the 2004 HESS data (resp. by the Egret data), when Egret (resp. HESS) 
data place a stronger constraint. Furthermore, we indicate for reference, again with a dashed black 
line, the exclusion limits which could have been drawn out of the 2003 HESS data only: this gives 
an idea of the improvement on the constraints on DM physics gained with the 2004 HESS data. All 
the exclusion limits have been computed in a model independent way, picking the most conservative 
option between a gauge boson dominated final state and a heavy quark final state. 

Obviously, the HESS data only constrain neutralino masses larger than the minimal energy probed 
by the experiment, respectively ~ 150 GeV for the 2004 data and ~ 280 GeV for the 2003 
data. The shape of the regions excluded by Egret and by HESS highlights a nice complementarity 
between satellite-borne experiments and ACTs in the exploration of the viable supersymmetric pa- 
rameter space. This complementarity will be further strengthened by the forthcoming space-borne 
GLAST and AMS experiments [89,90] and by the new generation ACTs [91]. We also remark that, on 
fully general grounds, current gamma rays data only place constraints on models providing a WMAP 
relic abundance (or on models with a rescaled DM density) if the DM density profile is cuspier than a 
NFW profile. Finally, since it is hard to imagine the DM density profile in the Galaxy to be shallower 
than the cored Burkert profile (see also Fig. the left panel shows that WIMP annihilation cross 
sections larger than 10^^^ cm'^s"^ are in any case in conflict with the Egret data. 

7 Conclusions 

For the ease of the reader, we collect below the main results of our analysis. 

• Adopting a fully numerical approach, we determined the range of WIMP masses and pair annihi- 
lation cross sections needed to fit the ACT data on gamma rays from the galactic center in terms 
of DM annihilations. The 2004 HESS data are not compatible with the hypothesis of a WIMP 
annihilating into a pure gauge bosons final state as well as that of a Kaluza-Klein DM candidate; 
the option of a mixed T'^r^-bb final state, viable within supersymmetric models, including some 
discussed in this paper (Sec. 13.21) . remains open, as a purely DM-annihilation interpretation of 
the HESS data. In the conservative "6esi spectral functions'^ approach, the WIMP annihilations 
interpretation of the HESS data is still valid, at 90% C.L., in the mass range 6 < m^^-ZTeV < 30 
(Sec© 

• We determined the maximal neutralino mass compatible with a thermal neutralino relic abun- 
dance within the WMAP range in the context of two benchmark minimal models, mSUGRA 
and mAMSB. The largest possible masses we find are not compatible with the mass range of 
the HESS data, while an interpretation of the Cangaroo-II data in the context of the mAMSB 
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model would require an extremely cuspy profile and an enhancement mechanism of the resulting 
thermal neutralino relic abundance (Sec. 13. Ij) 

• Enlarging the parameter space of mSUGRA and of niAMSB to general soft supersymmetry 
breaking Higgs mass terms, we showed that multi-TeV neutralinos are cosmologically allowed, 
and predicted to produce a detectable signal at future ton-sized direct detection experiments. 
Models with non-universal Higgs masses provide large LSP masses and pair annihilation cross 
sections, compatible with the ranges needed to fit the ACT data, and suitable gamma-rays 
spectral functions (Sec. 13. 2() 

• We gave an analytical and a numerical estimate of the largest possible neutralino pair annihilation 
cross section in the general minimal supersymmetric extension of the standard model, both for 
models with and without gaugino mass unification. The largest cross sections in the MSSM 
occur for resonantly annihilating neutralinos with maximal gaugino-higgsino mixing (Sec. |^ 

• The inclusion of non-perturbative electro-weak resonant effects can produce pair annihilation 
cross sections larger than the maximal bounds we derived, but only for model-dependent narrow 
ranges of the neutralino mass. The resulting gamma rays spectral function, moreover, is not 
compatible with a DM annihilation origin for the HESS data (Sec. 14. 1[) 

• We then addressed the issue of the largest possible neutralino mass in the MSSM, compatible with 
a "WMAP" thermal relic abundance. We studied the consequences of a full non-perturbative 
treatment of QCD effects in the pair annihilation cross section of a viable neutralino coannihilat- 
ing partner, the gluino. We showed that in this context, considering extreme non-perturbative 
scenarios, a neutralino with a mass much larger than 100 TeV, hence well beyond the s-wave 
unitarity limit, can yield a low enough thermal relic abundance (Sec. |SJ 

• We numerically assessed the maximal neutralino pair annihilation cross section and the largest 
possible "supersymmetric factor" ({av) /m'^), relevant for any indirect DM detection rate, for 
models with a WMAP relic abundance, or with a "rescaled" cross section (see Eq. ()31() ). The 
results we presented provide a conservative best-case scenario for indirect supersymmetric DM 
detection. As a particular instance, we showed how available data on the flux of gamma rays 
from the galactic center, as measured by Egret and HESS, constrain the viable supersymmetric 
parameter space (Sec. ^ 
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